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Abstract. We review our work on the p and a resonances derived from the Inverse Amplitude 
Method. In particular, we study the leading 1 /N c behavior of the resonances masses and widths and 
their evolution with changing m n . The 1 /N c expansion gives a clear definition of qq states, which 
is neatly satisfied by the p but not by the a, showing that its dominant component is not qq. The 
m K dependence of the resonance properties is relevant to connect with lattice studies. We show that 
our predictions compare well with some lattice results and we find that the pn% coupling constant 
is m % independent, in contrast with the an% coupling, that shows a strong m K dependence. 
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Light hadron spectroscopy lies beyond the realm of perturbative QCD. At low ener- 
gies, however, one can use the QCD low energy effective theory, named Chiral Perturba- 
tion Theory (ChPT) [1], to describe the dynamics of the lightest mesons. ChPT describes 
the interactions of the Goldstone bosons of the QCD chiral symmetry breaking, namely, 
the pions, by means of a effective lagrangian compatible with all QCD symmetries in- 
volving only the pion field. The infinite tower of terms in this lagrangian is organized 
as a low energy expansion in powers of p 2 /A 2 , where p stands either for derivatives, 
momenta or masses, and ~ 4nf n , where f n denotes the pion decay constant. ChPT 
is renormalized order by order by absorbing loop divergences in the renormalization of 
higher order parameters, known as low energy constants (LECs), that parametrize the 
high energy QCD dynamics and carry no energy or mass dependence. They depend on 
a regularization scale /i but after renormalization the observables are independent of this 
scale. The value of the LECs depend on the underlying QCD dynamics and are deter- 
mined from experiment. Up to the desired order, the ChPT expansion provides a sys- 
tematic and model independent description of how observables depend on some QCD 
parameters like the light quark mass m = (m u + m c i/2) or the number of colors, N c [2]. 

The use of ChPT is limited to low energies and masses, nevertheless, combined 
with dispersion relations and elastic unitarity it leads to a successful description of 
meson dynamics up to energies around 1 GeV, generating resonant states not originally 
present in the lagrangian, without any a priori assumption on their existence or nature. 
In particular, we find the p and a resonances as poles on the second Riemann sheet 
of %% elastic scattering amplitudes. With this approach we can then study some of 
these resonances properties, like their spectroscopic nature through their mass and width 
dependence on N c , or their dependence on the pion mass in order to connect with lattice 
studies. In the following sections we review this "unitarized ChPT" approach, named 



the Inverse Amplitude Method (IAM) [3, 4, 5], and then apply it to study the leading 
l/Nc behavior and the chiral extrapolation of the p and a mesons. 

The p and o resonances appear as poles on the second Riemann sheet of the (7,7) = 
(1,1) and (I, J) = (0,0) nn scattering partial waves of definite isospin, I and angular 
momentum J, respectively. Elastic unitarity implies for these partial waves, t(s), and 
physical values of s below inelastic thresholds, that 

Imt(s) = a(s)\t(s)\ 2 lml/t(s) = -o(s), with a(s) = 2p/ v G, (1) 

where s is the Mandelstam variable and p is the center of mass momentum. Conse- 
quently, the imaginary part of \/t is known exactly. However, ChPT amplitudes, being 
an expansion t ~ t 2 + U H , with tt = 0(p k ), can only satisfy Eq. (1) perturbatively 

lmt 2 (s)=0, lmt 4 (s) = a(s)t 2 (s) ... (2) 

The resonance region lies beyond the reach of standard ChPT. This region however, can 
be reached combining ChPT with dispersion theory through the IAM [3, 4, 5]. 

The analytic structure of the nn scattering amplitude t(s), consisting on a right cut 
extending from s t f, = 4m 2 to <*>, and a left cut from — <*> to 0, allows to write a dispersion 
relation for the auxiliary function G(s) = f|Cs)/f(.s) 

G(s) = G(0) + G , (0> + ±G"(0)s 2 + - f° ds' J mG(/) . , +LC(G)+PC, (3) 

L 7i J Sth s fi (s' — s — ie) 

where the integral over the left cut has been abbreviated as LC(G) and PC stands from 
the pole contributions in the scalar wave corresponding to the Adler zero. The terms in 
Eq.(3) are evaluated using unitarity and ChPT as follows: The right cut (RC) is exactly 
evaluated considering the elastic unitarity conditions Eqs. (1), (2): ImG(s') = —hates'). 
The subtraction constants only involve the amplitude and its derivatives evaluated at 
s = 0, so they can be safely approximated with ChPT: G(0) ~ t 2 (0)-t 4 (0), G'(0) ~ 
t 2 (0) - ^(O), G"(0) ~ -^(O). The LC, being suppressed by 1 /s' 3 (s' - s), is weighted at 
low energies, so it is appropriate to approximate it with ChPT: LC(G) ~ —LC{t 4 ). The 
PC counts 0(p 6 ), it has been calculated explicitly [6] and it is numerically negligible 
except near the Adler zero, away from the physical region. 

Neglecting PC for the moment, taking into account that t2(s) is just a first order 
polynomial in s, and that a dispersion relation can be also written for U, we can write 
Eq.(3) as 

G(5)~?2(0)-4(0>-?4(0)-4(0)5-i4(0> 2 -/?C(?4)-LC(? 4 ) = ?2^)-?4^), (4) 

which immediately leads to the IAM formula t IAM (s) = ^yz^t^ ■ The IAM formula 
satisfies exact elastic unitarity and, when reexpanded at low energies, reproduces the 
ChPT expansion up to the order used to approximate the subtraction constants and the 
left cut. Here we have presented an 0(p 4 ) IAM but it can be generalized to higher 
chiral orders. Note that in the IAM derivation ChPT has been always used at low 
energies, to evaluate parts of a dispersion relation whose elastic unitarity cut has been 
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FIGURE 1. Left: p and a l/N c scaling <9(//). Right: a l/N c scaling 0(/? 6 ). 

taken into account exactly. Thus, there are no additional model dependencies in the 
approach, which is reliable up to energies where inelasticities become important. Taking 
the pole contribution into account leads to a modified IAM formula [6] which is almost 
indistinguishable from the ordinary one except in the Adler zero region, where the 
modified formula should be used. Actually, we use the modified IAM in this work since, 
as it will be shown below, one amplitude pole gets near the Adler zero region. 

This simple IAM formula is able to reproduce KK scattering data up to roughly 1 GeV 
and generates the p and a poles with values of the LECs compatible with standard ChPT 
[5]. The 1 /N c expansion is implemented in ChPT through the LECs, whose leading 1 /N c 
scaling is known from QCD. Also, the m n dependence of IAM agrees with ChPT up to 
the order used. Hence, it is straightforward to study the leading l/Nc behavior and the 
m dependence of the resonances generated with the IAM. 

The QCD l/N c expansion [2] provides a clear definition of qq bound states: their 
masses and widths scale as 0(1) and 0{l/N c ) respectively. The QCD leading l/N c 
behavior of the ChPT parameters (f K , m n and the LECs) is well known. Hence, by 
scaling with N c the ChPT parameters in the IAM, the N c dependence of the p and a 
mesons mass and width has been determined [7, 8]. They are defined from the pole 
positions as ^/s^y = M — zT. Note that we should not take too large Nc values, since the 
N c — > °° is a weakly interacting limit, where the IAM approach is less reliable [9]. Also, 
for very large N c , a tiny admixture of qq in the physical state could become dominant, 
but this does not give any information about the physical state dominant component. 

Fig. 1 (left) shows the p and o mass and width N c scaling. It can be seen that the p 
follows remarkably well the expected behavior of a qq state, confirming that the method 
obtains the correct N c behavior of well known qq states. In contrast,the o does not follow 
that qq pattern, allowing us to conclude that its dominant component is not qq. 

Loop contributions play an important role in determining the o pole position. Since 
they are l/N c suppressed compared to tree level terms, it may happen that for larger 
N c they become comparable to tree level 0(p ) terms, which are subdominant in the 
ChPT series, but not N c suppressed. Thus we checked the 0(p 4 ) results with an 0(p 6 ) 
IAM calculation [8]. We defined a # 2 -like function to measure how close a resonance 
is from a qq behavior. First, we used it at 0(p 4 ) to show that it is not possible to find a 
set of LECs that makes the o to behave predominantly as a qq state. Next, we obtained 
an 0(p 6 ) data fit where the p qq behavior was imposed. Figure 1 (right) shows the M a 
and r CT N c scaling obtained from that fit. Note that both M a and T a grow near N c = 3, 
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FIGURE 2. Top Left: Movement of the a (dashed lines) and p (dotted lines) poles for increasing 
m n on the second sheet. The filled (open) boxes denote the a (p) pole positions at m % = 1, 2, and 
3 x m^ ys , respectively. Top Right: Comparison the IAM M p dependence on m K with some recent lattice 
results[14]. Bottom Left: Comparison of the p (light) and cr (dark) mass dependence on m n . Bottom 
Center: Comparison of the p (light) and cr (dark) width dependence on m K . The dotted (p) dot-dashed 
(a) lines show the decrease due to only phase space assuming a constant coupling to %%. Bottom Right: 
p and a couplings calculated from the pole residue. In all panels, the bands cover the LECs uncertainty. 



confirming the 0(p 4 ) result of a non qq dominant component. However, for N c between 
8 and 15, where we still trust the IAM, M a becomes constant and r a starts decreasing. 
This may hint to a subdominant qq component, arising as loops become suppressed as 
N c grows. Finally, by forcing the o to behave as a qq, we found that in the best case this 
subdominant component could become dominant around N c > 6 — 8, but always with an 
N c — > °° mass above 1 GeV instead of its physical ~ 450 MeV value. This supports the 
emerging picture of two low energy scalar nonets, one of exotic nature below 1 GeV and 
another of ordinary qq nature above 1 GeV. 

ChPT also provides an expansion of m K in terms of m (at leading order m\ ~ m). 
Thus, by changing m n in the amplitudes we see how the IAM poles depend on m. We 
report here our analysis of the p and a properties dependence on m n [10]. 

The values of m n considered should fall within the ChPT applicability range and allow 
for some elastic regime below KK, that would almost disappear if m n > 500, which 
would be the most optimistic applicability range. We expect higher order corrections to 
be more relevant as m n increases. Thus, our results become less reliable as m n grows. 

Fig. 2 (top left) shows the evolution of the a and p pole positions as m K is increased. In 
order to see the pole movements relative to the tzk threshold, which is also increasing, we 
use units of m K , so the threshold is fixed at y/s = 2. Both poles move closer to threshold 
and they approach the real axis. The p poles reach the real axis at the same time that they 
cross threshold. One of them jumps into the first sheet and becomes a bound state, while 
its conjugate partner remains on the second sheet practically at the very same position 



as that in the first. In contrast, the o poles go below threshold with a finite imaginary 
part before they meet in the real axis, still on the second sheet, becoming virtual states. 
As m K increases, one pole moves toward threshold and jumps through the branch point 
to the first sheet staying in the real axis below threshold. The other a pole moves down 
in energies away from threshold and remains on the second sheet. Similar movements 
were found within quark models [12] and a finite density analysis [13]. 

Fig. 2 (top right) shows our results for M p dependence on m % compared with some 
lattice results [14] and the M p PDG value. In view of the incompatibilities between 
different lattice collaborations, we find a qualitative good agreement with lattice results. 
The M p dependence on m n agrees also with estimations for the two first coefficients of 
its chiral expansion [15]. 

In Fig. 2 (bottom left) we compare the m n dependence of M p and M a , normalized to 
their physical values. The bands cover the LECs uncertainties. Both masses grow with 
m K , but M a grows faster than M p . Above 2.4m^ hys , we show two bands since the two a 
poles lie on the real axis with two different masses. 

In the bottom center panel of Fig. 2 we compare the m n dependence of T p and r,j 
normalized to their physical values: note that both widths become smaller. We compare 
this decrease with the expected phase space reduction as resonances approach the %% 
threshold. We find that Y p follows very well this expected behavior, which implies that 
the pnn coupling is almost m K independent. In contrast, r CT deviates from the phase 
space reduction expectation. This suggests a strong m n dependence of the o coupling 
to two pions, which we confirm with a explicit calculation of the resonances couplings 
from the pole residues as shown in the bottom left panel. 
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